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^ ■ Abstract 



We discuss the computational complexity of the perturbative evaluation of scatter- 



ed ' ing amplitudes, both by the Caravaglios-Moretti algorithm and by direct evalua- 

tion of the individual diagrams. For a self-interacting scalar theory, we determine 
^ . the complexity as a function of the number of external legs. We describe a method 

^ ! for obtaining the number of topologically inequivalent Feynman graphs contain- 

(^ I ing closed loops, and apply this to one- and two-loop amplitudes. We also com- 

^ ■ pute the number of graphs weighted by their symmetry factors, thus arriving at 

exact and asymptotic estimates for the average symmetry factor of diagrams. We 



present results for the asymptotic number of diagrams up to 10 loops, and prove 
that the average symmetry factor approaches unity as the number of external legs 



Qj ' becomes large. 



1 Introduction 

With the advent of high-energy colliders such as LHC and TESLA, high-multiplicity 
final states will become ever more relevant, increasing the need for efficient eval- 
uation of complicated multi-leg amplitudes. Performing such calculations by a 
direct evaluation of all relevant Feynman graphs is computationally hard in the 
sense that the number of graphs increases with N roughly as N !, the total number 
of external legs. For example, the 2^8 purely gluonic amplitude in QCD con- 
tains 10.5 million Feynman graphs at the tree level; and one may expect that loop 
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corrections (described by many more diagrams) will also be important. A compu- 
tational breakthrough has been achieved by the introduction of the Caravaglios- 
Moretti (CM) algorithm [ 1 1, in which the Schwinger-Dyson (SD) equations of the 
theory, rather than their decomposition in individual Feynman diagrams, are em- 
ployed, thus leading to a complexity of order c*^ (where c is a constant). Such 
methods, however, have to date only be formulated for the Bom approximation. 
Barring a revolutionary new method for solving the SD equations including loop 
effects'^, the most straightforward approach would seem to use the vertices of the 
effective, rather than those of the bare, action. In such an approach the effective 
vertices with up to N legs have to be employed, which increases the complexity 
of the CM algorithm. In order to assess the relative merit of the CM algorithm, 
it is therefore relevant to compare the computational complexity of the CM ap- 
proach to the number of higher-order Feynman graphs. In section 2, we calculate 
the number of individual diagrams, not weighted by their symmetry factors, in 
zero-, one- and two-loop level for four models of a self-interacting scalar theory. 
We also give the number of one particle irreducible graphs, that is needed in the 
sequence. In section 3 we give the number of diagrams, now weighted by their 
symmetry factors, for the four models, as they occur directly from the path in- 
tegral. Section 4 contains asymptotic estimates, in the number of external legs, 
for weighted and unweighted graphs. In section 5 we proceed in calculating the 
computational complexity of the CM algorithm in one and two loops. In section 6 
we compare the efficiency of the CM algorithm to that of the individual-diagram 
approach. 

2 Counting diagrams 

We consider a self- interacting scalar theory with arbitrary vertices of the type cp*^, 
k = 3, 4, We define the 'potential' 

V((p]=}2e^^ , (1) 



k>3 



k! 



where e^ is 1 if the <p^ interaction is present, otherwise it is zero. We shall spe- 
cialize to a number of cases: 

cp^ theory : V((p) = (p^/6 , 



'^In informal discussions, all the experts agree that this would be a tremendous advance — but 
no-one has a clue on how to approach it. 



(p theory 

gluonic QCD 

effective theory 



V((P) = (p724 , 

V((p) = (pV6 + cpV24 , 

V(cp)=e'P-1-(p-(pV2 ; (2) 



but ahernative theories are easily implemented. 

2.1 Counting tree diagrams 

Tree diagrams can be conveniently counted by means of the SD equation. This 
hinges on the fact that, at the tree level, all diagrams have unit symmetry factor. 
The counting of tree diagrams has been described in detail in 0|3]lll, and here we 
briefly recapitulate these results. Let us denote by a[n) the number of Feynman 
tree diagrams contributing to the 1 -^ n amplitude, and define the generating 
function 

Pictorially, we denote this as 



*o(x] = ^p (4) 




By considering the alternatives when entering the blob from the left, we easily see 
that 

*o(x) = Y_ £i<+i — ^-^ffl, , (5) 

k>2 ^^^ ' 

where the right-hand side contains k blobs. This implies that (^o[x) obeys the 
equation 

ct)o(x)=x + V'(ct)o(x)) . (6) 

Since V((p) is of order (9((p^), this SD equation can easily be iterated starting with 
(t»o(x) = 0, and the desired <x{yv) can be read off once the iteration has proceeded 
far enough. Notice that 

V"(4)o(x)) = 1--l- , v(p)(cl)o(x)) = -l--^V(^-^'(*o(x)) (p>3) , 

(7) 



so that the higher derivatives of V((|)o(x)) are completely expressed in terms of 
4)o(x) and its derivatives. 

The asymptotic behaviour of a(n) for large n is determined by the singularity 
structure of <^q[x]. Since <^q[x] cannot reach infinity for finite values of x, the 
singularities take the form of branch cuts, where c()o(x) remains continuous but 
(as it turns out in all cases studied so far) ^'q[x] diverges. We have 



x = cl)o-V'(ct)c 



dx 
dct)o 



l-V"(4)o) 



(8) 



and the dominant singularity is reached for that value cf)c for which V"((|)c) = 1 
and 

Xc = 4)c-V'(4)e) (9) 

is closest to the origin^. This value is always located on the positive real axis, 
where cj)o(x) is concave and monotically increasing for x < Xc. Taylor expansion 
then gives the structure of the branch cut: 



4)o(x) ~ 4)c 



2Xr 



from which we conclude that, for large n, 

, , / x;^ n! 

an 






(10) 



c 



n! 



27rV(3)(ct)c] n3/2x^ 



An yv^/HT"^ 



(11) 



with C = a/2/V(3'(4)c)- In the table we give the relevant numbers for the four 
case theories. 



theory 


4)c 


Xc 


C 


(p3 

gQCD 

effective 


1 

-1+V3 
log(2) 


1/2 

\/8/9 

V3-4/3 

21og(2]-1 


v/2 

21/4 

(4/3)1/4 
1 



^Here, we disregard the possibility that there are several such values, arising from a symmetry 
of the potential such as in the case of theories with only cp™^ interactions (m > 4). These cases 
are treated in detail in |6| and references therein. The asymptotic results given here are 'coarse- 
grained' . 



2.2 Counting one-loop diagrams 

When closed loops are introduced, an SD-type equation itself cannot be used to 
count the number of topologically inequivalent graphs. This stems from the fact 
that the SD-type equations are local in the sense that they only consider (in a re- 
cursive manner) what happens at a single vertex of a diagram, while the topology 
of a graph containing closed loops is a global property of the whole graph. In- 
stead, one has to settle for an order-by order and topology-by-topology treatment. 

Every one-loop diagram can be viewed as a single closed loop, to which tree- 
diagram pieces (which we call leaves) are attached. From 




k>l 



ek+2 , , ^k 



V"(ct)o(x)) 



(12) 



where the sum has k blobs again, and we have introduced the shorthand notation 
V, we see immediately that the number of one-loop graphs can be completely 
expressed in terms of v. The generating function of Li(rL), the number of all 
one-loop non-vacuum graphs with precisely n external legs, is given by attaching 
leaves to a closed loop in all possible ways: 



-Cilx) = Y — Li(rL 

^-- TV. 



rL>1 




(13) 



The standard combinatorics for collecting the various external legs into leaves, and 
inspecting the symmetry properties of the resulting graphs, show that a one-loop 
graph with ra leaves has precisely the 'natural' symmetry factor l/(2ra), with 
two important exceptions: the graphs with one or two leaves have an additional 
symmetry since, for the one-leave graph, the loop line may be flipped over, and for 
the two-leave graph the two internal loop lines may be interchanged. This leads 
us to the strategy for computing the number of topologically inequivalent graphs: 

• Write down the vacuum graphs, with their 'natural symmetry factor'; 



Attach leaves in all possible places; 



• Multiply by the order of the residual symmetry left over after the particular 
attachment. 

Performing this program for the one-loop case, we find 



CAx] 



1 1 
1 

1 



^+r'+L 



Tn>3 
1 



1 

2ra 



1 111 

-V + -V -^logll -v) 



^v + ^v^+^logl^oM) 



(14) 



The number of one-loop diagrams with n external legs is given below for some 
theories 



N 


cp^ 


CP^ 


gQCD 


effective 


1 


1 





1 


1 


2 


2 


1 


3 


3 


3 


7 





14 


15 


4 


39 


7 


99 


111 


5 


297 





947 


1,104 


6 


2,865 


145 


11,460 


13,836 


7 


33,435 





167,660 


209,340 


8 


457,695 


6475 


2,876,580 


3,711,672 


9 


7,187,985 





56,616,665 


75,461,808 


10 


127,356,705 


503,440 


1,257,154,920 


1,730,420,592 



2.3 Counting two-loop diagrams 

At the two-loop level, there are three topologically different vacuum diagrams. 
These are: 



a: 



OO^ • "O^ • '^oos • 



(15) 



where we have indicated their 'natural' symmetry factor. Since these graphs con- 
tain vertices, we must also accommodate leaves attaching themselves to vertices: 




v*^'(*W)-vt^'(0) 




V'4'((t)(x))-V<4'(0) . (16) 



In case no leave happens to be attached, the expression for the vertices read, of 
course, V'^'(O) and V*'*'(0), respectively. This prohibits, for instance, the occur- 
rence of a three-point vertex in a cp"* theory. For each of the graphs we have to 
admit the possibility of zero, one, or more leaves on each line, and that of leaves 
on any vertex. For the determination of the residual symmetries it must be re- 
membered that lines without leaves on them may be interchanged, and vertices 
without leaves may be interchanged, provided the 'anchoring' of the graph to the 
external legs contained in every leaf present permits such an interchange. As a 
simple example, the vacuum graphs themselves, without any leaves on them, have 
a residual symmetry of precisely 8,12, and 8, respectively, so that indeed they will 
be counted precisely one time. For graph a there are now 2 x 3^ = 1 8 cases to 
be considered, and for b and c we have 2^ x 3^ = 1 08 cases. The results for their 
generating functions are: 



^(a)i 



C\-\x] 



(b). 



C 



4'^(x) 



^v'^H*) 



1 

12 



1 +v 



1 



1 -V 



+ 4 



+ (v'3'(0))^f2 + 3(l+v) + (l+v)^ 



1-v)3 



V'^'(*))'^fl+v+ ^ 



1 -V 

2 



1 -V 



+4(V(3'(0)) (1 +v 



(17) 



The total number L2(rL] of two-loop graphs with precisely n external lines is 
therefore given via 



Cy{x] 



^-- n! 



C^^\x) + C^^\x) + C^^\x) . 



(18) 



TL>0 



Below, we give again the results for our specific theories. 



N 


cp^ 


CP^ 


gQCD 


effective 





2 


1 


3 


3 


1 


3 





6 


7 


2 


10 


3 


29 


35 


3 


58 





217 


273 


4 


465 


42 


2,214 


2,876 


5 


4,725 





28,365 


38,034 


6 


57,900 


1,485 


436,780 


604,320 


7 


829,080 





7,847,420 


11,202,156 


8 


13,570,515 


97,335 


161,048,720 


237,187,552 


9 


249,789,015 





3,715,400,500 


5,645,523,408 


10 


5,105,239,650 


10,210,200 


95,156,789,700 


149,180,360,320 



The extension to three or more loops is a matter of establishing the vacuum 
diagrams. For the three-loop case, however, there are 15 such graphs. Dressing 
them with leaves leads to a larger number of cases to be considered, ranging from 
54 to 1 1,664 per graph. 

2.4 Counting amputated diagrams 

Loop diagrams containing tadpoles or seagulls are constant contributions to lower 
order diagrams and are usually ignored. Moreover, diagrams containing self- 
energy loops on external legs are absorbed, during the renormalization procedure. 
Removing such diagrams from the above results is a simple task. One has to sub- 
stract all contributions from (a) diagrams with loops carrying zero or one vertex 
and (b) diagrams carrying two vertices one of which is connected with an external 
leg while the other is a single propagator. 

For the one-loop case one has to substract the first graph in eq. [T21 as a set of 
tadpole or seagul diagrams, as well as a contribution from graphs of the form 



<:^j 



With these modifications the generating function reads 

1 



/:^{x] = -jv + -v^ 



■log(1-v)-xcl)o(x)V>^'(0) 



(19) 



The number of amputated one-loop diagrams for our test theories is given 
below: 



N 


cp^ 


^>^ 


gQCD 


effective 


3 


1 





4 


4 


4 


12 


3 


39 


43 


5 


117 





437 


502 


6 


1,290 


75 


5,800 


6,916 


7 


16,425 





90,450 


111,660 


8 


239,400 


3675 


1,627,640 


2,077,944 


9 


3,944,745 





33,258,715 


43,883,696 


10 


72,627,030 


303240 


761,405,820 


1,037,955,824 



For two-loops diagrams one has to consider separately each vaccuum graph. 
All graphs containing loops with less than two vertices should be removed, as well 
as a variety of special cases which lead to non-amputated diagrams. 

The generating functions for each of the three basic topologies becomes 



4"V) 



c?\.) 



4'V) 



8 ^^'\ (l-v)2 



1 

U 



V(3'(*c 



-V<4)(0)V<^'(0)xc|)o(x) 
2/6-Mv + 9v^-2v^' 



~{y^'\o) 



Mx(^Q + 



1-V)3 

6 - 3v2 + 2v3 + v4' 



V(4'(0)) xct)o(x)-V'3'(0)V(^'(4)o)-V'^'(0)Vt4)(0)2x4)o(x) , 



V'^'l*) 



2v2(2-vl 



:i-v)3 



lvt3H*]v<3)(o)x^^t^ 

8 1 — vH 



+ VUO] - 



2 4 x^ 



81 -V 



(20) 



The exact number of two-loop connected amputated diagrams for our test the- 
ories is given below 



N 


^>^ 


CP^ 


gQCD 


effective 


3 


4 





28 


37 


4 


63 


9 


457 


600 


5 


870 





7,285 


9,760 


6 


12,945 


460 


128,675 


177,160 


7 


212,940 





2,552,165 


3,617,824 


8 


3,874,815 


35,315 


56,538,055 


82,588,784 


9 


77,605,290 





1,387,411,690 


2,089,438,256 


10 


1,700,078,625 


4,090,800 


37,407,699,175 


58,096,995,744 



2.5 Counting IPI diagrams 

The same methods as above can easily be employed in order to compute the num- 
ber of one-particle irreducible (IPX) diagrams. We simply restrict ourselves to 
the IPX vacuum bubbles; and, since IPX diagrams cannot have any vertex in their 
leaves, we simply replace 4)o(x) in the arguments of V", V'^', V'"^', ... by x. For 
the generating function of the IPX one-loop diagrams, we therefore have 



c]^'M 



1 



1 



-W+-W 



^logd 



■ w 



w 



v 



(21) 



The resulting numbers are given in the following table. 



N 


cp^ 


CP^ 


gQCD 


effective 


1 


1 





1 


1 


2 


1 


1 


2 


2 


3 


1 





4 


5 


4 


3 


3 


12 


17 


5 


12 





57 


83 


6 


60 


15 


390 


557 


7 


360 





3,195 


4,715 


8 


2,520 


315 


30,555 


47,357 


9 


20,160 





333,900 


545,963 


10 


181,440 


11,340 


4,105,080 


7,087,517 



At the two-loop level, we similarly find 



c?'M 



lv^'\x) 



1 +W + 



1 — w 



+ 4 



10 



yO) 



X 



+ 



1 -W)3 



12 ' ' [ 1 -w 



(22) 



Numbers are given below. 



N 


cp^ 


CP^ 


gQCD 


effective 





1 


1 


2 


2 


1 


1 





3 


4 


2 


2 


2 


9 


13 


3 


7 





40 


62 


4 


36 


12 


265 


410 


5 


240 





2,230 


3,499 


6 


1,860 


225 


22,485 


36,213 


7 


16,380 





261,135 


435,852 


8 


161,280 


8,295 


3,418,695 


5,944,000 


9 


1,753,920 





49,712,670 


90,309,029 


10 


20,865,600 


481,950 


794,102,400 


1,510,208,963 



3 Counting with symmetry factors 

The counting of diagrams including their symmetry factors is a somewhat simpler 
task, which can be performed on the basis of the path integral itself. In |5| this 
has been discussed in detail. However our approach here is somewhat different. 
One can expand the generating function of the number of connected diagrams per- 
turbatively around cp = and get a series in x (the source). Or, alternatively, one 
can expand perturbatively around the tree level one-point function cp = ^q. This 
shift eliminates the source x in favour of the tree level one-point function (|)o(x), 
and reveals the vaccuum graph dressing procedure that we employed above. 

3.1 Counting diagrams with symmetry factors 

Consider the generating function for the number of disconnected diagrams of a 
scalar theory with arbitrary couplings and a source x: 



Zfx] 



N 



dcp exp 



1,1 2 



V((p]+x(p] 



(23) 



11 



with N = 1 /v^ItS. Expanding around the tree level approximation (^o of the 
one-point function, i.e. setting cp ^ cj^o + cp, and making use of the Schwinger 
-Dyson equation for c|)o(x) gives 



Z(x)=N exp(-ls(ct)o] + ^) 
n, n. 



dcp e-^^t'P' (24) 



with 

S(^)=l^r(*£)^^_fv'"'(4.o)^ (25) 

2 ^ — n! 

n=3 

The generating function of the number of connected diagrams is then 



W(x) =?Llog(Z(x)] = -S((l)o) +xct)o+Tilog(N 



dcpe^^^t'P') (26) 



We see that it can be seen as a sum of the tree level part plus higher order correc- 
tions. These corrections can be written as the generating function for the vacuum 
diagrams of a theory with action S ( cp ) . The Feynman rules corresponding to this 
action can be read off directly : 

• J^vk^] = *o for every propagator. 

• V'^H't'o) for every n-point vertex. 

Given the potential V( (p) of the theory one can expand the vertex terms in the 
exponential of eq.|2B calculate the Gaussian integrals and arrive at an expression 
for W(x) that contains only V"{(^o[x)) and its derivatives. In this way, given 
the tree level one-point function of the theory, one finds the number of graphs 
weighted by their symmetry factors to arbitrary order. 

Writing W(x) in an H expansion 

W[x) = Wo(x) +nWi (x) +n^W2(x) + . . . (27) 

and, performing the integral and collecting together the terms of the same order 
in H, we see that the one loop diagrams are generated by 

We can also find the generating function for the two loop diagrams 

Wi ) = l_^^i!!i*o^ , 5 V'3)(4)o)V(3)(4)o) 

^^^W 8(l-V"(4)o])2^24 (1-V"(cj)o)P ^ ^ 

12 



The factor ^ in front of the first term is the symmetry factor of the only 2-loop 
vacuum diagram with a 4-vertex (see fig [TSla). The factor -^ = -^ + j2 is the sum 
of the symmetry factors of the two vacuum diagrams with two 3-vertices (see fig 
El.b and .c)^ 

Writing the derivatives V'^^H'I'o) in terms of derivatives of cj^o (which can be 
done by differentiating the Schwinger-Dyson equation for (^q) one arrives at 

1 C 1(cl)'^)2 



'^'^''^ 8 (4)'o)2 6 (4)'o)3 
Below, we give results for our four case theories in 1 loop. 



(30) 



N 


cp^ 


^>' 


gQCD 


effective 


1 


1/2 





1/2 


1/2 


2 


1 


1/2 


3/2 


3/2 


3 


4 





15/2 


8 


4 


24 


7/2 


57 


63 


5 


192 





1,149/2 


658 


6 


1,920 


80 


7,230 


8,568 


7 


23,040 





218,175/2 


133,676 


8 


322,560 


3,815 


1,919,190 


2,430,816 


9 


5,160,960 





77,146,125/2 


50,484,016 


10 


92,897,280 


31,0940 


871,927,770 


1,178,963,856 



The results for the four theories in 2 loops are again collected below. 



N 


cp^ 


CP^ 


gQCD 


effective 


1 


5/8 





31/24 


17/12 


2 


25/8 


2/3 


25/3 


19/2 


3 


175/8 





1,777/24 


527/6 


4 


1,575/8 


149/12 


5,057/6 


1,037 


5 


17,325/8 





280,735/24 


44,726/3 


6 


225,225/8 


1,535/3 


1,149,515/6 


252,734 


7 


3,378,375/8 





86,813,545/24 


14,808,232/3 


8 


57,432,375/8 


111,755/3 


464,096,885/6 


109,143,424 


9 


1,091,215,125/8 





44,344,732,495/24 


8,085,390,392/3 


10 


22,915,517,625/8 


12,672,800/3 


292,590,237,275/6 


73,514,104,288 



^In fact one could even avoid performing the integral since the generating function for N loops 
is simply the sum of the vacuum graphs with N loops weighted by their symmetry factors using 
the Feynman rules for the S{(p) action given above. However, this presupposes that one knows 
what the symmetry factor of the specific vacuum diagram is. 



13 



Both in the one and two loop cases an intriguing pattern of denominators is appar- 
ent for large N values, which seems to persist (we have checked this for N up to 
50). 

The above procedure can easily be extended to higher-loop amplitudes as well, 
but since we have not computed the unweighted diagram sums we defer this dis- 
cussion to the case of asymptotically large N . 

3.2 counting IPI graphs 

The generating function for the one particle irreducible diagrams of a theory 
weighted by their symmetry factors can be obtained by the same prescription by 
substituting cj^o = x. Now, however, we have to take into account only the IPI 
vacuum diagrams. In the one loop case the only vacuum graph is IPI and the 
generating function is 

In the two loop case we have to take into account the vacuum graph with one 4- 
vertex (see figure [TSla) and only one of the two vacuum graphs with three vertices 
(see figure \T5\h) since the other one (see figure \T5\c) is not IPI . This alters the 
symmetry factor from ^ to y^- We get then 



W2(X) 



1 

8" 



V(4)| 



X 



1 



+ 



1 V'3)(x)V(3)(x) 



(32) 



V"(x))2 12 (1 -V"(x))3 

We give below the number of irreducible diagrams weighted by their symmetry 
factors in the 1-loop case for the four test theories : 





cp^ 


^>' 


(p^ + (p4 


effective 


N=l 


1/2 





1/2 


1/2 


N=2 


1/2 


1/2 


1 


1 


N=3 


1 





5/2 


3 


N=4 


3 


3/2 


21/2 


13 


N=5 


12 





57 


75 


N=6 


60 


15 


390 


541 


N=7 


360 





3,195 


4,683 


N=8 


2,520 


315 


30,555 


47,293 


N=9 


20,160 





333,900 


545,853 


N=10 


181,440 


11,440 


4,105,080 


708,7261 
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We give below the number of irreducible diagrams weighted by their symme- 
try factors in the 2-loop case for the four test theories : 





cp^ 


CP^ 


(p^ + cp"^ 


effective 


N=l 


1/4 





2/3 


19/24 


N=2 


1 


5/12 


41/12 


101/24 


N=3 


5 





89/4 


691/24 


N=4 


30 


21/4 


709/4 


5765/24 


N=5 


210 





1,660 


56,659/24 


N=6 


1,680 


135 


17,865 


64,0421/24 


N=7 


15,120 





217,035 


8,178,931/24 


N=8 


151,200 


5,775 


2,936,745 


116,422,085/24 


N=9 


1,663,200 





43,787,520 


1,827,127,699/24 


N=10 


19,958,400 


368,550 


713,163,150 


31,336,832,741/24 



4 Asymptotic estimates 



It is fairly easy to estimate the number of diagrams, both with and without their 
symmetry factors, for asymptotically large N. As before, the asymptotic be- 
haviour of these numbers is governed by the analytic structure of their generating 
functions close to that singularity which is closest to the origin (that is, around 
X ~ Xc). There, we have 



4>o(^) ~ 2^^^^^ 



-1/2 



C=f2/V(3'(4)c 



1/2 



(33) 



where Xc, (^c and C again depend on the theory. Let us first concentrate on the 
one-loop diagrams. Since v = 1 — 1/(|)q(x) has a square-root branch cut at the 
singular point, log( 1 — v) is more singular than v or v^, and we have 



1 / C 
£i(x) ~ -log --^= 

2 V^V^ 



X, 



Cf 



(34) 



We conclude that, for one-loop diagrams, the average symmetry factor of a given 
diagram is asymptotically equal to 1. The number Ki (N) of graphs contained in 
the one-loop N -point amplitude is asymptotically given by 



K,(N)4 ' ^' 



4 (xc)N N 



(35) 
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To illustrate the convergence of the weighted number of graphs to the unweighted 
number, we give the ratio of the coefficients of x^ in C\^ [x) to those of C-\ [x] 
as a function of N below, for the pure cp^ theory. The other cases show a similar 
behaviour^, in which the asymptotic regime is approached as 1 /-\/N : this can also 
be easily checked from the exact form of £i (x] close to the singularity. 



1.5- 



The asymptotic results for the higher-loop amplitudes can be established by the 
following reasoning. The leading contribution from each leave-dressed vacuum 
diagram is given by that part that has the highest degree of divergence as x ^ Xc. 
From each line in the vacuum graph, this is a factor 1/(1 — v) = (^'q{x). Fur- 
thermore, from each k-point vertex in the vacuum graph the leading contribution 
comes from the limiting behaviour of V''^'(ct)o(x)). Now, it is easily seen that, as 

X 7 Xc, 

Vt3'(*o(x))~-^ ^VM(4)o(x))~0 , k>4 . (36) 

We conclude that the leading behaviour of the number of unweighted graphs is 
given by those vacuum graphs that contain only three-point vertices. To get the 
number of unweighted diagrams at the L-loop level, therefore, we first compute 



^For the pure (p'^ theory, this holds in the 'coarse-grained' approximation |6l. 
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the normalized path integral for the pure cp^ theory, using the usual perturbative 
interchange between expansion of the potential term and integration: 



N 



L 

n>0 



dcp exp - 

i6n)l 



M- 2 
2^ 



A 



rCp-^ 



(2n)!(3n)!(576)^ \[i 



^' 



(37) 



The sum of all connected vacuum diagrams with interactions is then given by 

W = log(Z] , (38) 

in the expansion of which the L-loop contribution (L > 2) is given by the term 
with 'A^^^^. In this expression, it suffices to replace A by 2/C^ and (x by 1 /(^q[x). 
The result is 



W 



Y_ ^L c'-^ 



X 



i3(l-L)/2 



(39) 



L>2 



The first coefficients wl are given below. 



L 


Wl 


L 


Wl 


2 
3 
4 
5 
6 


5/48 

5/64 

1105/9216 

565 / 2048 

82825 / 98304 


7 

8 

9 

10 

11 


19675/6144 

1282031525/88080384 

80727925/1048576 

1683480621875 / 3623878656 

13209845125/4194304 



The asymptotic result for Kl(N), the number of unweighted diagrams contribut- 
ing to the L-loop n-point amplitude is therefore given by 



kl(n: 



r(N + |(L 



1 



xPc)'"V(|(L 



llxN 



(40) 



For the number of L-loop graphs weighted by their symmetry factors we may 
employ the following formulation of the SD equation: 



{np,q}>0 



1 



1 



T^P.q 



p,q>0 ^^P-qJ 



:q + i 






(41) 



17 



where the bracketed superscripts denote derivatives, and 

^ (p + q)np,q = L , m=1+_^(q + l)np,q . (42) 



p,q p,q 



The successive expressions for 4)l(x) in terms of lower-loop ones can straight- 
forwardly be worked out. For L = 1 , 2 these have been given in the previous 
section. If we now put in the approximate form of (^q{x) given in Eq. dTOb . it is 
easily checked (at least up to L = 10) that expression W is reproduced. Note 
that in this approximation the fourth and higher derivatives of V((|)o] vanish, so 
that Eq. lHTT) is actually more complicated than need be: nevertheless, by using the 
next-to-leading expression 

ct)o(x)~4)c-C(xc-x]^/2-C'(xc-x) , (43) 

it can also be checked that, indeed, the subleading behaviour of ^o{x) shows up 
only in the subleading terms in Kl(N ). We conclude that as N —) oo, the average 
symmetry factor of any Feynman diagram approaches unity. 



5 Complexity of the Caravaglios-Moretti algorithm 

5.1 introduction 

The CM algorithm, as first explicitly given in lHJ (and earlier implied by 0), 
consists of the computation of subamplitudes with one off-shell leg, the other legs 
corresponding to on-shell external legs of the transition matrix element. For tree 
diagrams, these subamplitudes can be unambiguously specified by the particu- 
lar set of external momenta involved because of momentum conservation. For 
detailed descriptions, we refer to lH] 0: here, we are only interested in the 
combinatorics of the algorithm. 

5.2 Complexity for tree level computations in any theory 

We assume an N-particle process, and set K = N — 1. Each subamplitude can 
then be encoded by a binary string with N bits, each referring to a given external 
particle. The bit is set to 1 if its external leg is involved in the subamplitude, and 
to otherwise. For instance, the string (1 , 1 , 0, 1 , 1 , 0, 0, 0, . . . , 0, 0) denotes that 
subamplitude in which the external particles with labels 1,2,4 and 5 are combined. 
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(0,0,1,1,1, 




using the vertices of the theory, into a single off-shell momentum. By the same 
convention, a string with a single 1 refers to the Feynman rule for a single external 
particle (a spinor or antispinor for fermions, a polarization vector for vector parti- 
cles, etctera). The CM algorithm combines subamplitudes into successively more 
complicated ones, culminating in the string (1 , 1 , 1 , . . . , 1 , 1 , 1 , 0), which, after 
multiplying with the external factor (0,0,0, ...,0,0,0,1) gives the final answer 
for the amplitude. It is clear that of the N external particles, one can be left out 
of the combinatorics since it has to be included only at the very end. The combi- 
natorial problem is, therefore, to determine the number of ways to decompose a 
string of K bits. An example of a possible decomposition is 

(1,0,0,0,0,0) 

(1,1,0,0,0,0» (0,1,0,0,0,0) 

(1,1,1,1,1,0) ^^ (0,0,1,0,0,0) 

(0,0,0,1,0,0) 
(0,0,0,0,1,0) 

where we have indicated the strings corresponding with the external legs and 
the various subamplitudes. The possible decompositions depend on the theory 
in question: the presence of an (ra + 1 ) -point vertex in the theory allows for a de- 
composition into ra smaller strings. In this paper, we shall only deal with theories 
of a single self- interacting field (gluonic QCD being an example): extensions to 
more fields are fairly straightforward. In recent implementations such as HELAC 
([|9i|), this decomposition can be recognised explicitly. 

Let us first consider a subamplitude's string with n I's being decomposed 
into ra smaller strings, each with at least one 1. This happens when, in the SD 
equation, an (ra + 1 ) -point vertex is encountered. The number of inequivalent 
decompositions, denoted by Crri(TT^), is given by 

Cnx(n) = -^ y — ^ , (44) 

where, of course, n-\ + Uj + ■ ■ ■ + rim = n. Note that the above equation as- 
sumes that all the subamplitudes containing ni,rL2, . . . ,rLm, external momenta 
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exist. This is always the case when a cp^ interaction is present in the theory^ but 
it is not true for a pure cp"^ theory for example. Then one has to introduce a factor 
that cancels the terms coming from combinations of non-allowed subamplitudes. 
We, nevertheless , proceed with our program to find a generating function for 
effective theories that always contain a three point vertex. We find 

Now, out of all bit strings of size K, there are precisely K!/rL!(K — n)\ strings 
containing precisely TL I's. The total number of decompositions involving (ra+l)- 
point vertices is therefore 

fm(K) = }1 ( M c^(n) , (46) 

n>0 ^ ^ 



SO that 



g^[x)^y^i^[K) = ^e^{e^-^r . (47) 

^ — K m 



K>o K! m! 



In the simple case of a pure cp^ theory we therefore have 

Q2[x]^Ue'^- le^^ + e'') = }^ ^ 1 (3X - 2'<+^ + 1 ) , (48) 

K>0 ■ 

SO that the number of decompositions necessary to arrive at an N-point amplitude 
is given by 

6 2^2' 

For a theory with both cp^ and cp"^ interactions such as gluonic QCD, we find a 
total of 

J_4N _ I2N , 1 

24 4 3 

decompositions. In QCD at the tree level, an improvement is possible. We can 
decompose the gluonic 4-vertex into two 3-vertices by employing an auxiliary 
field, as explained for instance in [6|. This brings the complexity down from 4*^ 

'^because then there is always the possibiHty of constructing a subamplitude containing Uk 
external momenta by combining a subamplitude containing n\c_ — 1 momenta with an external 
momentum in a three point vertex 
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to 3 , a worthwile improvement for large N. It is not to be expected, however, 
that this will be possible in higher orders. The effective action, therefore, will 
contain (ra + 1 ) -vertices for all ra > 2, and the generating function is therefore 



F(x) = 2_ QmM = exp (e" — 1 + x) — exp(2x) 



(49) 



Tri>2 



Below we give the number of decompositions, 

D(N) = ^gUN-l) , 



(50) 



Tri>2 



for not-too-large values of N . 



N 


D(N) 


N 


D(N) 


3 


1 


8 


4,012 


4 


7 


9 


20,891 


5 


36 


10 


115,460 


6 


171 


11 


677,550 


7 


813 


12 


4,211,549 



For asymptotically large values of N , we have to study the analytic structure of 
F(x). Since this function is analytic for finite x, D(N) must increase with N 
slower than N!. On the other hand, D(N) increases faster than c*^ for any finite 
c, which is reasonable since as N grows larger and larger values of ra come into 
play. This is also evident from the fact that the standard Borel transform of the 
series Ffx], 



dy e ^ T{xy] 



1 



+ 



1 -2x 
does not converge for any positive value of x. 



dij exp (— -y + e"^ — 1 + xy] 



(51) 
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20 



40 



60 



80 



100 



The above plot shows the behaviour of the ratio (log(N)/N)(D(N)/D(N - 1)] 
as a function of N for 3 < N < 1 00. For high N, this ratio is approximately (but 
not quite) a constant. 

5.3 Complexity in one and two loops 

Consider a general theory with ra-point vertices. Each subamplitude of level n 
(containing n specific external momenta), can be constructed by combining two 
or more lower-level subamplitudes in a three- or more-point vertex. When using 
an (ra+ 1 )-point vertex the subamplitude is built by ra lower level subamplitudes 
and the number of different ways for this to happen is given by eq.l44l 

Each term in the series represents the number of ways to construct the sub- 
amplitude of level n using subamplitudes of level ni , . . . , nra- The computational 
cost of each such subamplitude involves, (assuming that there is an ra+ 1 vertex in 
the theory) contributions from the following posibilities: All lower subamplitudes 
are free of loop corrections and the vertex is an ordinary one (this gives the tree 
level subamplitude)^. It can also be that one of the subamplitudes contains already 

^That is provided that the lower level subamplitudes exist! This always happens when the 
theory involves cp^ interactions. In the pure cp^ theory, however, we have to modify the calculation 
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a loop correction (that occured in previous steps in the CM. algorithm) and the 
vertex is an ordinary one (see fig 15^ . The subamplitudes containing loop correc- 
tions can, however, be of level 2 or higher since the level one subamplitudes are 
the external legs which we consider amputated. There are, therefore, tu— ^|6i_tvi 
different possibilities. Finally there is the case that all subamplitudes are free of 
loop corrections but the vertex is actually a loop (see the last term in fig 15^ . The 
number of different possibilities is now equal to the number of IPI diagrams with 
one loop and ra + 1 legs, which we denote with Jm,,!- 




+ 




+ 



+ 




(52) 



The cost of computing the specific subamplitude via an ra + 1 vertex is there- 



fore 



L 



n! 



ra! ■* — rii In?! . . . riTT, ! 

TLl ,...,nm ' ^ "^ 



6(^n,-n)(V^(l+m-_^5i,,J + J^,i) (53) 



where we have included a factor Vm, = 1 if the ra + 1 vertex is in the theory 
and Vra = if not. The cost of the subamplitude is then found by summing over 
ra. There are ^|Lj.,_ ,, different subamplitudes. The computational cost of the 
whole algorithm in units of effective vertices is then 




(VU1+m-^5i,^] + U (54) 



where K = N — 1 . In the following table we present the results for the four test 
theories. 



to exclude combinations where one of the n^'s is equal to 2 since in such a theory there are no 
level 2 subamphtudes. 
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^>^ 


^>' 


cp^ + cp" 


effective 


N=l 














N=2 














N=3 


2 





5 


6 


N=4 


18 


4 


46 


57 


N=5 


114 





340 


442 


N=6 


720 


105 


2,715 


3,713 


N=7 


5,368 





26,346 


37,411 


N=8 


49,686 


3,395 


315,035 


459,056 


N=9 


553,766 





4,474,868 


6,688,320 


N=10 


7,112,700 


149,140 


72,741,355 


112,139,709 



In order to include the two loop correction one has to add to the above formula 
a term (ra — ^ 5i ^i ) (ttl — ^ Sm ~ 1 ) for the possibility that two of the lower 
subamplitudes have a one loop correction and a term equal to ra — ^ 61^4 for 
the possibility that one of the subamplitudes has a two loop correction. There is 
also the possibility that one of the lower subamplitudes is of 1-loop order and the 
vertex itself is a 1-loop IPI graph. This costs an extra term of Jm.,i (tu — ^ 5i tt^ ) 
Moreover one has to add the number of IPI graphs with 2 loops and ra + 1 legs, 
Jrn,,2- Hence we now have, writting S^ = Hi ^i.nt 



K 



n=2 ^ ' -n^l ni,...,nn 



n! 



ni ! . . . tItt^! 



-(A + B + C) 



where K = N — 1 and 

A = Vtt^(1 +m-Sni + (m-SnJ(Tn-Sni -1) +m-Sni; 

B = Jm,1 + Jm,l(TTL-SnJ 
C = Jm,2 

The results for the four test cases are presented below: 
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^>^ 


^>' 


cp^ + cp^ 


effective 


N=l 














N=2 














N=3 


9 





45 


68 


N=4 


102 


16 


566 


857 


N=5 


957 





6,414 


9,837 


N=6 


9,740 


610 


81,560 


127,451 


N=7 


114,677 





1,201,556 


1,920,824 


N=8 


1,546,986 


32,151 


20,211,345 


33,181,094 


N=9 


23,395,461 





380,938,056 


644,468,452 


N=10 


390,310,512 


2,574,670 


7,929,937,496 


13,861,514,611 



One should be aware of the fact that the above results are obtained under the 
assumption that the computational cost for every effective vertex that might in- 
clude one or two loop IPI graphs is the same. 

6 Comparison of the complexity of the CM. algo- 
rithm to the diagrammatic approach 

We present below the ratio of the computational complexity of the CM. algorithm 
over the number of diagrams one has to calculate in the customary diagrammatic 
approach, for our four test theories^". For each case the ratio for a calculation in 
tree, tree plus 1-loop, and tree plus 1- and 2- loop level is presented. 



complexity of CM. algorithm / number of diagrams 




cp^ 


^>' 




Lo 


Lo + Li 


Lo + Li + L2 


Lo 


Lo + Li 


Lo + Li + L2 


N=3 


1.00 


1.000 


1.500 


- 


- 


- 


N=4 


2.00 


1.200 


1.307 


1.000 


1.000 


1.231 


N=5 


1.666 


0.864 


0.955 


- 


- 


- 


N=6 


0.857 


0.516 


0.679 


2.00 


1.235 


1.119 


N=7 


0.318 


0.309 


0.498 


- 


- 


- 


N=8 


0.093 


0.199 


0.375 


1.575 


0.858 


0.819 


N=9 


0.022 


0.136 


0.286 


- 


- 


- 


N=10 


0.005 


0.095 


0.220 


0.636 


0.468 


0.584 



"only amputated, tadpole/seagull-free diagrams are considered 
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complexity of CM. algorithm / number of diagrams 




(p-^ + (p4 


effective 




Lo 


Lo + Li 


Lo + Li + L2 


Lo 


u + u 


Lo + Li + L2 


N=3 


1.000 


1.000 


1.363 


1.000 


1.200 


1.619 


N=4 


1.500 


1.070 


1.132 


1.500 


1.212 


1.325 


N=5 


1.00 


0.736 


0.828 


0.961 


0.837 


0.956 


N=6 


0.409 


0.451 


0.605 


0.381 


0.519 


0.691 


N=7 


0.121 


0.283 


0.454 


0.109 


0.327 


0.515 


N=8 


0.028 


0.190 


0.347 


0.025 


0.217 


0.392 


N=9 


0.005 


0.132 


0.268 


0.005 


0.150 


0.302 


N=10 


0.001 


0.094 


0.208 


0.001 


0.107 


0.234 



One should note that the CM. algorithm will actually perform better than de- 
picted by the above numbers, when compared with the straightforward diagram- 
matic approach, since we consider the cost of a step in the CM. algorithm (i.e. the 
calculation of a subamplitude which corresponds to the calculation of an effective 
vertex) equal to the cost of the computation of a whole diagram. That is the reason 
for the apparently poor performance of the CM. algorithm in the case of tree level 
cp"* theory. 

We, therefore, conclude that the Caravaglios-Moretti algorithm is more effec- 
tive than the straightforward diagrammatic approach, in the tree as well as the one 
and two loop level, by a factor that increases rapidly with the number of external 
legs, even though this increase is less rapid in the one- and two- loop level than in 
tree level. 
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